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This paper describes a new measurement of the flux ratio of positive and negative muons from cosmic-ray 
interactions in the atmosphere, using data collected by the CMS detector at ground level and in the under- 
ground experimental cavern. The excellent performance of the CMS detector allowed detection of muons in the 
momentum range from 3 GeV to 1 TeV. For muon momenta below 100 GeV the flux ratio is measured to be 
a constant 1.2766 ± 0.0032(stat) ± 0.0032(syst), the most precise measurement to date. At higher momenta 
an increase in the charge asymmetry is observed, in agreement with models of muon production in cosmic-ray 
showers and compatible with previous measurements by deep-underground experiments. 
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1. Introduction 

This paper describes a measurement of the ratio 
the flux of positive to negative muons produced by 
cosmic ray interactions in the atmosphere, using the 
CMS detector [l|, Q . Cosmic ray particles (predomi- 
nantly protons) interact with the upper atmosphere, 
initiating a hadronic shower. Kaons and pions pro- 
duced in these showers may decay to muons, which 
can then reach the Earth's surface. These muons are 
therefore an easily detectable signal of a cosmic ray 
interaction, and their properties provide information 
on the initial cosmic ray particle and development of 
the hadronic shower in the atmosphere. 

The momentum distribution of muons has been pa- 
rameterised Q based on the kinematics of primary 
cosmic-ray particles and of decays of secondary par- 
ticles, and this parametrisation can be used as the 
base to extract the muon charge ratio, R, as described 
in 1: 
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where the constants f n and fx are the fractions of 
all pions and kaons that are positive, and must be 
determined from data. The constant r\ sets the rel- 
ative pion and kaon fractions in cosmic-ray showers 
and the accepted value of 0.054 [3j is taken. There is 
also a dependence on the vertical component of the 
muon energy, E^cos0 z , where 9 Z is the zenith angle. 
Finally, there are two energy scales: 115 GeV and 
850 GeV, above which pions and kaons respectively 
have a greater probability to interact rather than de- 
cay in the atmosphere. Between 115 and 850 GeV 
therefore, the importance of kaons as a source of 
muons is expected to rise. As the value of fx is larger 
than /„■ due to the preference of K + A production over 
K + K~~ production, the muon charge ratio is also ex- 
pected to rise in the energy range between approxi- 
mately 100 GeV and 1 TeV. Previous measurements 
of the muon charge ratio 0, S S 0, H S E3, EH have 




Figure 1: Two cosmic ray muons reconstructed by CMS. 
Left: using only the muon detectors; right: using the 
muon and tracking detectors. 



been sensitive to energies below or above this expected 
rise, and this paper describes the first measurement to 
bridge this gap, while also being the most precise to 
date in the range 5-100 GeV. Results are presented in 
terms of both momentum and the vertical component 
of the momentum, p.cos0 z . 

The CMS detector [l2|, one of the experiments at 
the CERN Large Hadron Collider 13], is designed 
to study proton-proton collisions at energies up to 
14 TeV, and search for new physics at this energy 
frontier. The detector is described in detail else- 
where [14|, with the main feature being a supercon- 
ducting solenoid magnet with a diameter of 6 m, pro- 
viding a field of 3.8 T. Inside the magnet are a central 
tracking detector, surrounding the luminous region, 
and a crystal electromagnetic calorimeter and brass- 
scintillator hadronic calorimeter. Muons are measured 
in gas-ionisation detectors embedded in the steel re- 
turn yoke for the magnet. The ability to accurately 
measure TeV-energy muons is central to the physics 
goals of the LHC programme, with the added benefit 
that CMS also makes an excellent high-energy cos- 
mic muon detector. Figure [T] shows two typical cos- 
mic muon events reconstructed by CMS: using only 
the muon detectors, and a combination of muon and 
tracking detectors. 
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The data used in this analysis were collected 
in two phases of the commissioning of CMS: the 
2006 "Magnet Test and Cosmic Challenge" (MTCC) 
runs and the 2008 "Cosmic Runs at 4 Tesla" 

(CRAFT08) pH. These two datasets were analysed 
separately due to the different detector conditions, 
and the results then combined. Much of the work in 
designing the analyses focused on ensuring the cancel- 
lation of many charge-dependent systematic and ac- 
ceptance effects, simplifying the analysis and reducing 
the overall uncertainties. This was complicated by the 
fact that, unlike some other experiments, CMS does 
not change the polarity of the magnet. 



2. MTCC Analysis 

During the MTCC runs, a part of the CMS de- 
tector was installed at ground level in the assembly 
hall above the main experimental cavern, at 509 m 
above sea level, 46°18.57' north, 6°4.62' east. While 
two sections of the muon system were available, for 
this analysis a fiducial volume was denned using only 
one of these sections, to ensuring a symmetric accep- 
tance for positive and negative muons, as shown in 
Fig. [2] A total of 330,000 muons were selected for the 
analysis, based on a trigger using the MB2 and MB4 
layers of the detector (see Fig. |5J). As the MTCC 
analysis was at the Earth's surface, it was possible to 
detect muons with momenta down to 3 GeV. Recon- 
structed muons were corrected for energy lost travers- 
ing the detector, and the charge ratio as a function 
of p and p.cos0 z calculated. Finally, a further correc- 
tion for charge mis-identification is derived from sim- 
ulation (using the CORSIKA Monte Carlo [17] and 
GEANT [l8[ detector simulation) to obtain the final 
distribution. 



3. CRAFT08 Analyses 

For the CRAFT08 runs, the full CMS detector 
was installed in the experimental cavern, 89 m below 
the surface. To identify cosmic ray muons traversing 
the detector (rather than muons produced in proton- 
proton collisions at the centre of the detector), spe- 
cial triggers and reconstruction algorithms were devel- 
oped 0, • Two analyses were based on the data 
collected in these runs: one using "global" muons, and 
one using "stand-alone" muons. Global muons require 
the muons pass through the CMS muon system and 
central tracking system, which greatly improves the 
momentum resolution while greatly reducing the ac- 
ceptance, relative stand-alone muons which use the 
muon system only. A total of 300,000 events were se- 
lected for the global analysis, and 3 million for the 
stand-alone analysis. In both cases, selections on the 
number of hits and the fit quality of the muons are 
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Figure 2: Fiducial region for the MTCC analysis, 
showing the charge-symmetric acceptance and trigger for 
muons. The solid red and blue lines show muons that 
would be accepted, the dashed lines rejected. 



used to minimise charge confusion, and events with 
more than one reconstructed muon rejected. Figure [3] 
shows the charge ratio reconstructed in the detector 
for both analyses. 

Both analyses then follow a similar method to ex- 
tract the final result: the reconstructed muon tracks 
were individually propagated to the Earth's surface, 
correcting for energy loss both in the detector and in 
the rock above the cavern (the average energy loss be- 
ing around 50 GeV). However, leading down to the 
cavern are three access shafts, as shown in FiglU and 
even very low energy muons can reach the detector by 
traveling down these shafts. As two of the shafts are 
off the central axis of the detector, they break the ac- 
ceptance symmetry for positive and negative muons: 
coming from one side of the detector, the solenoid 
magnet will bend positive muons into the detector ac- 
ceptance, and negative muons out. To overcome this, 
these shafts were excluded from the analysis, along 
with symmetric regions on the other side of the detec- 
tor to preserve symmetry. 

The distributions of p and p.cos9 z at the surface 
are then formed, and corrected for detector resolution 
and charge confusion effects. The detector resolution 
for global muons was determined in data, using "split 
tracks" . Muon tracks crossing the whole detector were 
split (see Fig. [T]), and the resolution determined from 
the difference in the measurements of the incoming 
and outgoing legs. For stand-alone muons, a GEANT 
based detector simulation was used to determine the 
resolution, which was cross-checked by comparing the 
stand-alone to central tracker resolution in the global 
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Figure 3: The reconstructed charge ratios from the two 
CRAFT08 analyses, as a function of transverse 
momentum at the centre of the detector (pt CA ) 
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Figure 4: Layout of the CMS experimental cavern, 
showing two of the three access shafts. 



muon data sample. 



4. Systematics 

For the MTCC analysis, the uncertainty on the de- 
tector alignment is the main limiting factor, at around 
2%, increasing with muon momentum to around 
10%. Much smaller uncertainties due to charge mis- 
identification and knowledge of the magnetic field also 
contribute. For the CRAFT08 analyses, systematics 
were assessed due to the uncertainty on the material 
above the CMS cavern and any charge-dependent ef- 
fect on the muon rate, the impact of the event selec- 
tion and trigger, and knowledge of the magnetic field. 



However, the largest single source of uncertainty was 
again the understanding of the detector alignment, 
which rises to 3-4% at high momentum. 

The detector was aligned using cosmic muons dur- 
ing the MTCC and CRAFT08 runs however as 
well as uncertainties on this alignment, there remains 
the possibility of "weak modes" : deformations of the 
detector which do not change the fitted track % 2 , but 
do affect the measured track momenta. As these weak 
modes are \ 2 invariant, the standard alignment pro- 
cedure is to some level blind to them, and the main 
concern for this analysis is the class of weak modes 
which induce a (signed) offset in the measured curva- 
ture, therefore shifting the measured momentum for 
positive and negatively charged muons in opposite di- 
rections. To constrain such weak modes, the "end- 
point test" was developed. The number of measured 
muons should fall to zero as the measured momentum 
goes to infinity, or as q/pr goes to zero. Fitting the 
measured q/pr distribution can test for any shift of 
this minimum away from zero, and set limits on the 
size of any weak modes. In the CRAFT08 analysis, 
an offset of 0.043 ± 0.022 TeV^ 1 was found, the domi- 
nant uncertainty at high momentum. It is worth not- 
ing that such weak modes will be further constrained 
by LHC collision data, as more alignment techniques 
become possible. 



5. Combination 

The fully corrected MTCC, global and stand- 
alone analyses were combined using a standard pro- 
cedure [221 ] . Previous measurements have shown 
that the charge ratio is a constant in the range 5- 
100 GeV, so in a first combination all data points be- 
low p.cos8 z — 100 GeV are combined into a single 
measurement. This involves 6 points from the MTCC 
and 3 points from each of the global and stand-alone 
analyses and a full treatment of correlations between 
data points and analyses, yielding: 

R = 1.2766 ± 0.0032(stat) ± 0.0032(syst) (2) 

This is the most accurate measurement to date of the 
charge ratio in the range 5-100 GeV. 

Then, the data points from different analyses in 
each bin were combined, again keeping track of all 
correlations. The resulting combination is shown in 
Fig. [5J and the charge ratio as a function of p.cos6 z 
is also compared to previous measurements. The 
CMS results agree well with the L3+C data @ be- 
low 400 Ge^and with the Utah [5], MINOS 
and OPERA (ll| measurements above 400 GeV. The 
CMS data are also consistent with a rise in the charge 
ratio above 100 GeV, and fitting with Eq. [1] yields 
f n = 0.553 ± 0.005 and f K = 0.66 ± 0.06, with a 
fit x 2 /d-o.f. = 7.8/7. It is worth noting that while 
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Figure 5: Top: the three measurements of the charge 
ratio, R, and combination as a function of momentum, 
with the global and stand-alone analysis points offset for 
clarity. Bottom: the CMS combination compared to 
previous measurements as a function of p.cos9 z . 



the L3+C data contribute many points in the range 
20-400 GeV, these points are in fact limited by a cor- 
related systematic uncertainty. As a result, the CMS 
measurement is also the tightest constraint to date on 
/tt, which dominates the value of Eq.[T]below 100 GeV. 

6. Conclusions 

CMS has made a new measurement of the charge 
ratio of atmospheric muons. The measurement is the 



most precise to date in the range 5-100 GeV, and the 
first to span the interesting range 100-1000 GeV. It is 
worth noting that the dominant systematic uncertain- 
ties, arising from the knowledge of detector alignment, 
will be reduced by LHC proton-proton collision data, 
making future cosmic ray muon measurements from 
CMS an interesting prospect. 



References 



[i 

[2 

w 

[s: 

[7; 

[8 

[9 
[10 

[11 

[12 
[13 
[14 

[15 

[16 
[17 

[18 

[19 
[20 
[21 
[22 



CMS Collaboration, Phys. lett. B 692 (2010) 

83-104. 

CMS Collaboration, CMS Pvlg rMUO-lO-OOll 
(2010). 

T. Gaisser, "Cosmic Rays and Particle Physics" , 

Cambridge University Press 1990. 

P. A. Schreiner, J. Reichenbacher, and M. C. 

Goodman, Astropart. Phys. 32, Issue 1 (2009) 

61. 

G. K. Ashley, J. W. Keuffel, and M. O. Larson, 

Phys. Rev. D 12 (1975), no. 1 20. 

J. M. Baxendale, C. J. Hume, and M. G. 

Thompson, J. Phys. G: Nucl. Phys. 1 (1975) 

781. 

B. C. Rastin, J. Phys. G: Nucl. Phys. 10 (1984) 
1629. 

T. Hebbeker and C. Timmermans, Astropart. 
Phys. 18 (2002) 107. 

L3 Collaboration, Phys. lett. B 598 (2004) 15. 
MINOS Collaboration, Phys. Rev. D 76 (2007) 
052003. 

OPERA Collaboration, Eur. Phys. J. C 67 
(2010) 25. 

CMS Collaboration, JINST 3 (2008) S08004. 
J. Evans and P. Bryant, JINST 3 (2008) S08001. 
CMS Collaboration, CERN/IHCC 2006-001 

(2006) . ____ 
CMS Collaboration, CMS Note 2007/005 

(2007) . 

CMS Collaboration, JINST 5 (2010) T03001. 
D. Heck et. al., Forschungszentrum Karlsruhe, 
Report FZKA 6019 (1998). 
S. Agostinelli et. al, Nucl. Inst. Meth. A 506 
(2003) 250. 

CMS Collaboration, JINST 5 (2010) T03022. 
CMS Collaboration, JINST 5 (2010) T03002. 
CMS Collaboration, JINST 5 (2010) T03009. 
L. Lyons, D. Gibaut, and P. Clifford, "How to 
combine correlated estimates of a single physical 
quantity", Nucl. Inst. Meth. A 270 (1988) 110; 
A. Valassi, "Combining correlated measurements 
of several different physical quantities" , Nucl. 
Inst. Meth. A 500 (2003) 391. 



C30 



